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Role of deformation on giant resonances 
within the QRPA approach and the Gogny force 
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Fully consistent axially-symmetric-defornied Quasi-particle Random Phase Approximation 
(QRPA) calculations have been performed, in which the same Gogny DIS effective force has been 
used for both the Hartree-Fock-Bogolyubov mean field and the QRPA approaches. Giant resonances 

*^fi OC-"- OO oA *^n Oft 

calculated in deformed Si and Mg nuclei as well as in the spherical Si and Mg isotopes 

are presented. Theoretical results for isovector-dipole and isoscalar monopole, quadrupole, and 
octupole responses are presented and the impact of the intrinsic nuclear deformation is discussed. 

PACS numbers: 21.60.Jz, 24.30.Cz, 27.30.+t 



I. INTRODUCTION 

A major challenge in theoretical nuclear physics is the 
development of a universal approach able to describe 
the excited states of all nuclear systems with the same 
accuracy. From this perspective, models based on the 
Random-Phase- Approximation (RPA) [l| are well-suited 
in rigid nuclei as they describe on the same footing both 
individual and collective excited states. Such approaches 
have been successfully applied to different nuclear sys- 
tems. However, most of these first calculations were re- 
stricted to spherical nuclei with no pairing correlations. 
The generalization of these RPA calculations to all nu- 
clear systems requires to treat explicitly pairing corre- 
lations, the intrinsic nuclear deformation and the con- 
tinuum coupling. In practice, it is hard to achieve all 
these improvements at once. Choices have been made ac- 
cording to the physics under study. Pairing correlations 
have been included to describe excited states of open- 
shell nuclei, and then theQuasi-particle RPA (QRPA) 
formalism has been used ^ 0] . Because of the increase 
of the number of unstable nuclei experimentally accessi- 
ble, some effort has been made also to treat the con- 
tinuum coupling in RPA approaches 0, [^, Q as well 
as in QRPA ones [3, @, |3, UOj Ell- Progress has also 
been made in the development of high order (Q)RPA 
calculations [13, [ij, ^§\ ■ Since the work of Gering and 
Heiss 115 1 , d eformed (Q)RPA models have been devel- 
oped ia)ll3)lla]- Let us note that many of these improve- 
ments have also been applied on top of relativistic mean 
field calculations [1^, [20j [21|. It is worth pointing out 
that all these developments should be made along a con- 
sistent line: it has been shown that neglecting parts of the 
residual interaction in the RPA equations can strongly 
affect predictions [H, [H, [H [H, [ll, [13]. Then, more 
and more calculations are now performed using fully self- 
consistent (Q)RPA approaches [H, [2i|. 
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In this paper, a fully consistent axially-symmetric- 
deformed Hartree-Fock-Bogolyubov (HFB)-I-QRPA ap- 
proach using a unique effective force, namely the DIS 
Gogny interaction [29|, is presented. This approach can 
be applied in open or closed-shell nuclei and in spher- 
ical or deformed ones. The present study is the first 
application of this method and is devoted to giant res- 
onances in the s-d shell nuclei, ^^"^^Mg and ^®"^°Si. 
Many experimental results are available for light nuclei 
[13, [M [si, [H, [H, [H, [H [13, [H, [H, [40| and a recent 
reanalysis of experimental data in ^^Si shows that giant 
resonances in these nuclei are still of great interest [40] . 
These Magnesium and Silicon isotopes display different 
ground state deformations, then they are good candi- 
dates for studying the impact of intrinsic deformation 
on giant resonances. In the past, the shift in energy of 
isoscalar response for deformed nuclei was estimated with 
macroscopic models, as reviewed in section 4.7 of ref |41| . 

In section [Hi an outline of the QRPA method is pre- 
sented. Practical aspects of our approach are discussed 
in section IIIIl Mean field calculations for 22-28]y[g ^^^^ 
26-30gj g^j,g presented in section HVl together with results 
for intrinsic deformations and pairing energies. Isovector 
dipole responses are discussed in section |V| and isoscalar 
monopole quadrupole and octupole results in section FVll 
Results are compared with experimental data and sys- 
tematic laws. Conclusions are drawn in section fVIII 



II. FORMALISM 

Approaches based on the RPA have been found to be 
successful in explaining low-lying multipole vibrations 
and giant resonances. The QRPA extension, which al- 
lows one to treat on the same footing particle-hole (ph), 
particle-particle (pp) and hole-hole (hh) excitations is 
well suited for both closed-shell and open-shell nuclei. 
In particular QRPA approaches on top of HFB calcula- 
tions describe collective vibrations of the nuclear mean 
field and pairing field together with induced ground state 
correlations. 



In the QRPA approach, the quasi-boson operators 6*+ 
representing nuclear excitations are defined as: 



e. 






Xi^vM - Y^VjVi) , 



(1) 



where rj^ and rji are quasi-particle (qp) creation and an- 
nihilation operators, respectively. Xn and Yn are the 
amplitudes of the two quasi-particles (2-qp) excitations. 
The QRPA ground state |0 > and the QRPA excited 
states \n > satisfy the equations: 



6»+|0>=|n>, 6'„|0>=0. 



(2) 



QRPA equations are obtained through the second 
derivative of the energy functional built with an effective 
interaction, in the present work the Gogny force DIS [23 ]. 
QRPA equations can be converted into a matrix equa- 
tion: 



A B 
B A 



Xn 
Yr,. 



Xn 
-Yn 



(3) 



where a;„ are the energies of the QRPA excited states 
\n >. In Eq. ([3]), the matrices A and B are real and 
symmetric since the time reversal symmetry T and the 
TIi2 symmetry, where 112 is the reflection with respect 
to the (xOz) plane, are assumed. A and B are built using 
the same Gogny DIS effective force in all ph, pp and hh 
channels. 

To obtain the QRPA energies uon and the amplitudes 
Xn and K„, Eq. Q is rewritten as: 

{A -B){A + B) {Xn + Yn) = ujI {Xn + Yn) ■ (4) 

Introducing the matrix S which diagonalizes {A — B) 
with D = S^^{A-B)S, Eq. (g]) can be rewritten 

_pi/2_5.-i (^ ^ Q^ SD^'^D-^/'^S-^ {Xn + Yn) 

= UJ-'D-^/^S-'{Xn+Yn), (5) 

i. e. as the eigenvalue equation: 

CVn = LolVn, (6) 

where C and Vn are defined as: 

C ^D^'^S-^{A + B)SD^'^, (7) 

and 

Vn = D-'^^S-' {Xn + Yn) . (8) 

Let us note that since the matrices A and B are real 
and symmetric, S can be chosen orthogonal {S^^ = S ). 
Consequently C is a symmetric matrix. 

The QRPA energies w„ are then calculated as the pos- 
itive square-root of the eigenvalues oj^ of C. Using the 
same procedure for 

{A + B){A- B) {Xn - Yn) = Lol {Xn ~ Yn) , 



similar to Eq. ([4]), and the orthonormalization relation 
Xn Xn' — Yn Yn' = Snn' the QRPA amplitudes Xn and 
Yn are then obtained as: 



Xn 



i"« - 2 



i?l/2 



i^l/2 



D-'/'V^n 



-D-^/^V^n 



Vn 
Vn- 



(9) 



Such a procedure, used to solve the QRPA equations 
(Eq. ([3])), has the advantage of reducing the size of the 
matrix to be diagonalized by half. The condition of 
positive-definiteness of at least one of the matrices A + B, 
A — B can be relaxed using a generalized Cholesky de- 
composition [4a |. 

In the present approach, the matrix equations are 
solved in the finite space of N 2-qp excitations. Then, 
{A — B) and C are NxN matrices and are diagonalized 
using standard routines giving all eigenvectors and eigen- 
values. Let us note that uj^, the solutions of Eq. ([5]), can 
be negative and then lead to imaginary u;„ values. Such 
complex eigenvalues of the QRPA matrix represent un- 
stable collective modes. However, as discussed in Ref. [Ij, 
it is always possible to find a solution of the HFB equa- 
tions which makes these collective modes stable and con- 
sequently ensures that they have real frequencies. 

In the present work, axially-symmetric- 

deformed (ASD) HFB calculations are performed 
in even-even nuclei imposing some symmetries, T, 
TH2 and also axial and left-right symmetry. Then 
the projection K of the angular momentum on the 
symmetry axis and the parity tt are good quantum 
numbers. Consequently ASD-QRPA calculations can be 
performed separately in each K'^ block. In the following, 
the ASD-QRPA states \n) of Eq.© wiU be labeled with 
the angular momentum projection K and rewritten as 
\nK). Time- reversed states will be written \nK). 

States \JMK)n of good angular momentum J are ob- 
tained projecting the intrinsic excitations \nK) according 
to: 



\JMK)r 



^/2J+\ 



Att 



I dn (dIjk* {n) R {n) \nK) 



\J-i< nJ 



+ {-)■'-'' Di,_j, {n)R{n)\nK)\, (10) 



where J7 are Euler angles, Dj^^j^ Wigner rotation matrix 
elements, and R{^) the three dimensional rotation op- 
erator. 

In the QRPA approach all quantities are calculated rel- 
ative to their values in the ground state. The response 
of the system to an external field, such as a multipole 
operator Qx^, is obtained from |( QQQ\Q\fj_\JMK)\^ , i.e. 
the square of the transition matrix elements between pro- 
jected \JMK) states of Eq.HO]) and the projected QRPA 
ground state |000). Using Eq. (fTO|) . one gets: 

{{)m\Qx^j_\JMK)n - Yl ^JK{mju\nK)5Ku , (11) 
i/,|iy|<j 



where |0) is the ASD-QRPA vacuum of Eq.® and the 
Tjx coefficients are obtained from straightforward alge- 
braic calculations. 

In an axially-symmetric-deformed nuclear system, the 
response function of a given 3'^ contains different K'^ ~ 
C^, ztf^, ...± J'^ components. In spherical nuclei, all these 
components are degenerated in energy, then the response 
functions associated to any multipolarity can be obtained 
from K"^ = 0* results only. This property has been used 
to check the newly-built ASD-QRPA computer code. 

The collectivity of each QRPA state is ob- 
tained through its contribution to the total 
Energy- Weighted-Sum-Rule (EWSR) Mi{Qxf,) 
J2n=i^n\{000\Qx,.\JMK)^\\ The mean energy of the 
QRPA resonances in a given energy interval [E„ 
is calculated as: 



^nitn 7 ^max\ 



<E> 



[£„;„,£„ 



Ml [E, 



En 



rmn ; ^^max 



Mo[E, 



mm ■} ^max\ 



(12) 



where the moments M\ are 



Ma [E„,,n,Ernax] = J^ ^n\{000\Q x^^W M K) „\^ (13) 
n 

with Em.irr < UJ^r < ^r7 



-'min — *^n _ ^-'max- 



III. CALCULATION ASPECTS 

Before presenting the results on Magnesium and Sili- 
con isotopes obtained with the above formalism, let us 
discuss the practical aspects of the present calculations. 

First, in the present study, the HFB equations are 
solved in a finite Harmonic Oscillator (HO) basis. Con- 
sequently the positive energy continuum is discretized. 

Second, since, as discussed by Terasaki et al 26] QRPA 
consistency is much affected by inaccurate HFB results, 
the HFB equations have been solved with a very high 
degree of accuracy and extremely well converged quasi- 
particle states are used in QRPA calculations. 

Third, all HFB quasi-particle states are used to gen- 
erate the 2-qp excitation set (Eq.([T|)). This means that 
no cut in energy or in occupation probabilities is intro- 
duced. Our model space, includes 9 HO major shells. 
This is a large enough space for nuclei studied here. We 
checked that the calculated EWSR is altered by less than 
2% when going from 9 to 11 HO major shells. 

Fourth, as already emphasized, the omission in the 
QRPA calculations of parts of the two-body force which 
are included in the HFB mean field and pairing field may 
lead to large errors in strength functions. Such a lack 
of consistency is often responsible for a redistribution of 
the strength [H H^Jp, [gl, HE [g^] and a strong spuri- 
ous state mixing [2J]. In the present calculations, the 
same effective Gogny DIS nucleon-nucleon force is used 
to solve both the HFB and the QRPA equations in all 
ph, pp and hh channels. However, since the Coulomb 



exchange field is not taken into account in HFB, the cor- 
responding QRPA terms have been set to zero in order 
to maintain consistency. 

Finally, as RPA-type approaches restore broken sym- 
metries, violations of conservation laws by the HFB wave 
function lead to spurious states. Spurious states are ex- 
pected to have zero energy eigenvalues in the (Q)RPA 
spectrum [l| . In the present axially-symmetric-deformed 
QRPA calculations there are seven spurious states. Ex- 
cept for two of them related to particle number violation, 
spurious states differ from one another since no two of 
them have the same K'^ quantum numbers: K'^ = 0~, 
±1~ for translational ones, K"^ — ±1+ for rotational 
ones, and I-C" = 0+ for particle number ones. Since 
our calculations have been performed along a consistent 
line, with i) the same interaction in HFB and QRPA 
calculations, ii) highly converged HFB calculations. Hi) 
no cut-off in single particle excitations, and iv) a large 
enough HO basis, these spurious states are easily iden- 
tified with large values of A ~ —Y and they are found 
to be well separated from physical excitations. They are 
found at 0.001 keV energy for K'' = 0+, around 500 
keV energy for K'^ = ±1+ and less than 1.8 MeV for 
K^ ~ 0~,±1~. The spurious modes relative to transla- 
tional invariance are more sensitive to the finite size of the 
HO basis, which prevents us from fully restoring trans- 
lational invariance. In order to put all these spurious 
modes at zero energy, renormalization factors could have 
been introduced. However, as studied in Ref. [17], renor- 
malization factors associated to different spurious modes 
would differ from each other and no universal choice ex- 
ists. A corrected dipole operator could also be used but, 
as mentioned in Ref. l2^ responses for the corrected and 
uncorrected operators are indistinguishable in the case 
of fully consistent HFB-I-QRPA calculations. In view of 
this, the results presented below have been obtained with 
uncorrected operators and no renormalization factors has 
been included. 



IV. HFB RESULTS FOR 



^Mg AND 



■'Si 



Constrained ASD-HFB calculations have been per- 
formed in ^^^•^°Si and 22-28]y[g ^^^lei. Potential energy 
curves for the seven nuclei are plotted in Fig. [T] as func- 
tions of the axial deformation parameter (3: 



(3 



5 Ai?2' 



(14) 



where (720 is the mean value of the axial quadrupole op- 
erator Q20 = Y.i=i Vl67r/5 r^{i)Y2o{i), and i?^ = 

3/5(1.2^1/3)2 f^2^ 

We observe that i) the N = 16 nuclei, ^°Si and ^^Mg, 
are found spherical, ii) the Z = 12 nuclei, •^^Mg and 
2'^Mg, are strongly prolate with (3 = 0.47 and 0.52, re- 
spectively and Hi) Z and/or N= 14 ones, ^^Si, ^^Si and 
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FIG. 1: HFB potential energy curves V as functions of the axial deformation parameter /3 for 



J Si and 



'Mg. 



^^Mg, present large oblate deformations with (3 ~ -0.3. 
In all these nuclei, except spherical ones, the potential en- 
ergy curves exhibit a secondary minimum as function of 
axial deformation. HFB calculations in triaxial symme- 
try show that these minima are connected to the ground 
state through triaxial deformations. In these light nuclei, 
only two protons (or two neutrons) radically change the 
deformation of the HFB ground state: this is of course a 
consequence of shell effects, which are also responsible for 
the vanishing of the pairing energy in specific deforma- 
tion ranges. In open shell systems the pairing energy take 
on a large value. This value is smaller when there is a 
significant single particle gap above the Fermi level. The 
cancelling of the pairing energy is associated with closed 
shell systems 43]. Here the pairing energy is defined as 
Ep = -^Tr (Ak) and is calculated separately for neutrons 
and protons, with A the pairing field and k the abnor- 
mal density [43] . In Fig. [2l proton and neutron pairing 
energies are plotted as functions of /3 in dashed and solid 
lines, respectively. These pairing energies range from 
to 6 MeV. The most striking feature of Fig. [2] is that a 
strong competition is found between the Z, N = 12 pro- 
late sub-shells and the Z, N = 14 oblate ones. In ^2Mgi2 



and f4Sii4 nuclei, proton and neutron pairing energies are 
canceling at the same deformation: prolate for ^^Mg and 
oblate for ^®Si. On the contrary in ^^Si and ^^Mg proton 
and neutron pairing energies do not vanish at the same 
deformation. Corresponding potential energy curves of 
^^Si and ^^Mg (Fig. [T]) present two minima, one for each 
sub-shell. These minima are close in energy, in contrast 
to those in ^^Si and ^^Mg which are separated by more 
than 2 MeV. In the N = 16 ^"Si and ^^Mg nuclei, the 
neutron pairing energy almost vanishes for (3 ~ 0.3, and 
for (3 = and spherical shapes are stabilized. The N = 10 
neutron pairing in ^■^Mg presents a minimum at very large 
deformation. 

In the following, QRPA calculations are performed us- 
ing the quasi-particle states associated with the HFB so- 
lutions with minimum energy in all nuclei. These so- 
lutions are either prolate (^^~^^Mg) or oblate (^^~^^Si, 
^^Mg) or spherical (^^Mg, 3°Si). In ^'^Mg and ^^Si, since 
the total pairing energy is zero in the HFB ground state, 
QRPA calculations reduce to conventional RPA ones, 
where elementary excitations are of particle-hole type. 
First results for these two nuclei have been presented 
previously in a conference proceeding J44| . In the other 
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FIG. 2: Neutron (full lines) and proton (dashed lines) pairing energies Ep as functions of axial deformation parameter /3 in 
26-30Si and 22-28Mg. 



nuclei under study, the full QRPA formalism is neces- 
sary as either proton or neutron pairing correlations are 
present in the HFB ground state. 



V. DIPOLE RESPONSE 

As mentioned in [4l|, the first evidence of giant- 
resonance phenomena in nuclei was obtained in 1937 
[45| . and the notion of a dipole oscillation of the nu- 
cleus was recognized by Migdal in 1944 [46|. The re- 
sults of many years of experimental work on the Isovec- 
tor Giant Dipole Resonance (IVGDR) have been sum- 
marized and reviewed in several papers [43, |4^, |49J : 
in light nuclei the strength distribution is fragmented 
into several components, whereas it is separated into 
two major components in deformed nuclei. QRPA re- 
sponses to the isovector electric dipole operator Qio = 
J2f=i ^^^w in) - Y.'iLi ''^^lo in) are given in Fig. [3] as 
fractions of the EWSR. The K '^ = 0", iKr = 1" com- 



ponents are plotted in black and green, respectively. In 
the two spherical nuclei ^^Mg and '^°Si, only the K^ = 0^ 
component is drawn as |K \^ = 1~ states are degenerated 
with K^ = 0~ ones. In all nuclei, the predicted strength 
is spread over a wide energy range up to 50 MeV, such a 
broad width being consistent with the fact that these nu- 
clei lie in the middle of the s-d shell. Negligible strength 
is found at low energy which is in agreement with the fact 
that pygmy states are observed only in very exotic nuclei. 
In spherical nuclei ^^Mg and "^"Si, responses are similar: 
a highly collective state is found at 23 MeV and two less 
collective components around 20 MeV and 25 MeV. In 
deformed nuclei, the strength is found fragmented and 
splits up into mainly two components, as expected. In 
this study, the two components correspond to two dif- 
ferent angular momentum projections K. More precisely, 
in the prolate ^^"^''Mg nuclei the lower energy peak is 
associated to K = states whereas the higher one cor- 
responds to |K| = J = 1 states. The same features are 
found in the oblate nuclei ^^^•^"Si and ^^Mg, with lower 
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FIG. 3; (Color online) Fractions of the isovector dipole EWSR in ^^ ^®Mg and ^^ ^"Si. In deformed nuclei, K" = and 
|K 1'^ = 1~ components are indicated in black and green, respectively. 



energy parts composed of |K| = J = 1 states and higher 
ones composed of K = states. The energy split between 
K = and K = ±1 states is related to the intrinsic de- 
formation of the ground state with an opposite hierarchy 
for prolate and oblate shapes. In prolate nuclei the sym- 
metry axis is the long axis and the |K| = J component 
corresponds to the alignment of the angular momentum 
along this axis, whereas for oblate nuclei the symmetry 
axis is the shorter one and the K = component corre- 
sponds to the alignment of the angular momentum along 
the long axis. Such a splitting of giant resonance into two 
groups of transitions, corresponding to vibrations along 
the major and minor axis of a deformed spheroid was al- 
ready predicted in 1958 by Okamoto [5(| and Danos [51[ 
in a long range correlation model. Such an effect was 
used in order to explain part of the broadening of the 
resonance widths in non spherical nuclei. 

From another theoretical point of view, Bassichis at 
al. [54I have explained the GDR splitting in ^^Mg as 
due to configuration splitting: the high energy region 



(22 - 27 MeV) is formed by lp-> (ld-2s) transitions, 
whereas the lower one (13 - 22 MeV) contains mainly 
(ld-2s) -^ (lf-2p) transitions. Our results are compati- 
ble with this analysis: major components of the K^ = 0~ 
QRPA state at 19.35 MeV are ld5/2 -^ lf7/2 transitions 
with small contributions of 2p3/2 -^ 2sl/2 transitions, 
whereas major components of the [Kj'^ = 1~ QRPA 
state at 27.44 MeV are lp3/2^2sl/2 transitions with 
non-negligible components from ld5/2^2f5/2 ones. 

By comparing our results with experimental data [38| , 
a relatively good agreement is found in ^^^^^Mg. In 
Ref. dH the general structure of the ^^Mg IVGDR is de- 
tailed: a first group of resonances at 18 - 21 MeV forms 
the main peak of the GDR and a second one is a broad 
plateau at 21 - 27 MeV with structures at 23 MeV and 
25 MeV. The width of the GDR is found to be 9 ± 1 MeV, 
in fair agreement with our results. 

As mentioned in Ref. [41|, a systematic law can 
be used to estimate the mean energy of the IVGDR, 
E (IVGDR) = 31.2 A-^ + 20.6 A- 6 MeV. Systematic 



TABLE I: Theoretical (present work) and systematic mean 
energy values of the IVDGR (in MeV). 





E{GDR) MeV 






th. 


Syst. 


22Mg 


23.09 


23.4 


^^Mg 


23.18 


22.9 


^''Mg 


23.13 


22.5 


28Mg 


23.03 


22.1 


26Si 


22.91 


22.5 


28Si 


22.80 


22.1 


^«Si 


23.17 


21.7 



TABLE IL Theoretical (present work), systematic and exper- 
imental mean energy values of the GMR (in MeV) . 





E{GMR) 


MeV 






th. 


Syst. 


Exp. 


22Mg 


21.23 


28.5 




^^Mg 


21.06 


27.7 


21.0±0.6 


26Mg 


21.99 


27.0 




2SMg 


21.33 


26.3 




2«Si 


21.62 


27.0 




2«Si 


21.76 


26.3 


21.25±0.38 


«°Si 


21.85 


25.7 





values are compared to our results in Table HI Theo- 
retical values of the IVGDR are obtained by integrat- 
ing the strength from 15 to 35 MeV. They are found 
roughly constant irrespective of the mass number, con- 
trary to the values obtained from the systematics which 
vary from 21.7 MeV (^OSi) to 23.4 MeV (^^Mg). In our 
calculations, the resonance energy is found almost inde- 
pendent of the deformation of the nucleus, as already 
predicted by Danos j5i]. In all seven nuclei the IVGDR 
is predicted to exhaust around 85% of the EWSR. In pro- 
late ^^"^''Mg nuclei, high energy components dominated 
by |K|'^ = 1~ states are found to exhaust 54% of the 
EWSR. Low energy components, dominated by K = 
states, exhaust 29% of the EWSR. In the three oblate 
nuclei, high energy components, dominated by K = 
states, exhaust around 35% of the EWSR and low en- 
ergy ones around 50%. Let us note that in a systematic 
study [5^ Goriely and Khan have shown that an equally 
distributed strength between low energy and high energy 
components give optimal agreement with experimental 
El transitions. 



VI. ISOSCALAR MULTIPOLE RESPONSES 



A. Monopole response 



18 MeV and 25 MeV, the high energy one being very 
fragmented. The low energy component is found to ex- 
haust 40% (22Mg) and 35% (^^Mg) of the EWSR, and 
the high energy broad component 44% (^^Mg) and 48% 
(2%Ig). 

Third, only a few discrete levels emerge in spherical 
^*Mg and "^"Si nuclei: the low energy strength is concen- 
trated into two major peaks around 18 MeV and 22 MeV 
and a non negligible high energy component is found at 
30 MeV. 

In Table [Til theoretical predictions of resonance 
mean energies are compared with available data [36J, 
ISTJ . Theoretical values are calculated in the energy 
range [9 - 41 MeV] as was done for experimental 
data [Sa [331 . Values obtained from the systematic law 
E {GMR) = 80 A-i MeV are given in the same TablelH 
Comparing our predictions with systematics and experi- 
mental data we find that: i) the mean energy is here pre- 
dicted almost constant, irrespective of the mass of the 
nucleus, contrary to the A~^/^ dependence of the sys- 
tematic law, but ii) theoretical energies are smaller than 
systematics, QRPA results being close to experimental 
data in 2**Mg and ^^Si. 

The splitting of monopole responses for ASD nuclei is 
related to the coupling with K = Q component of the 
quadrupole responses ^41.] . 



Fig. m displays monopole responses of ^^~^^Mg, and 
26-30gj g^g fj-actions of the EWSR. The most striking fea- 
ture is that the seven nuclei under study present features 
that can be separated in three classes depending on the 
deformation of their ground states. 

First, the responses of the three oblate nuclei ^^Mg and 
26-28gj display a prominent peak around 20 MeV and a 
broad component at a higher energy exhausting less than 
30% of the EWSR. The low energy part of the resonance 
exhausts 65%, 62% and 66% of the EWSR in ^SMg, ^631 
and ^^Si, respectively. The experimental response of the 
oblate ^®Si [3a] displaying a sharp peak close to 18 MeV 
and a large fragmentation up to 35 MeV, is qualitatively 
reproduced by our results. 

Second, the responses of the two prolate nuclei, namely 
^^Mg and ^''Mg, display a two-peaked structure around 



B. Quadrupole response 

Since the strength in light nuclei is found to be highly 
fragmented, the experimental identification of various gi- 
ant resonances is difficult. As mentioned in references 
[20, [sill [33 J this led in the 70 's to substantial disagree- 
ments between various studies [33, [34, .35]. QRPA re- 
sults for quadrupole responses in the seven nuclei under 
study are displayed in Fig. [5l where K'^ = 0+, ±1+, ±2+ 
components are in black, green and purple, respectively. 
Discrete spectra have been folded with a Lorentzian dis- 
tribution, C {x) = r/27r {x^ + r2/4) with a F = 2 MeV 
width, and the result plotted as blue lines. 

In all seven nuclei, quadrupole responses are predicted 
to lie between a few MeV and 30 MeV without any 
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FIG. 4: Fractions of the monopole EWSR in 22-28]yjg ^^^^ 26-30gi_ 



strength above, contrary to dipole and monopole re- 
sponses which extend up to 50 MeV (see sections |V] and 
I VI Ap . Discrete low- lying states are found up to 18 MeV 
and a Giant Quadrupolc Resonance (GQR) is predicted 
around 20 MeV. Let us note that the first low-lying state 
is found to be |K|'^ = 2+ in all deformed nuclei except 
^^Si. The occurrence of |K|'^ = 2+ states at a very 
low energy is an indication of the softness of these s-d 
shell nuclei with respect to the triaxial degree of free- 
dom as mentioned in section IIVI As to monopole and 
dipole modes, in the oblate nuclei (^^~^*Si and ^^Mg) 
|K| = J = 2 components of the GQR have lower energies 
than K = ones, while opposite behaviors are observed 
in prolate nuclei (^^~^'*Mg). The energies of the K"^ = 1+ 
components are around 20 MeV in the two prolate Mag- 
nesium isotopes and close to 22 MeV in the oblate nuclei, 
that is in between the energies of the K'^ = 0+ and of the 
|K|'^ = 2+ components, in all deformed nuclei. The same 
hierarchy was found by Hagino 'Yq et al in ^"^Mg and in 
■^^Mg. As the splitting of the quadrupole resonance due 
to the deformation is not very large, only one broad peak 



is found in the GQR. The width of the GQR is found to 
be smaller than the ones of the giant dipole and monopole 
resonances. However the energy splitting between K = 
and |K| = 2 components is arround 5 MeV in agreement 
with the 12-14 A~^/'^MeV splitting obtained from macro- 
scopic models (see section 4.7 of ref (41|]). The quadrupole 
responses of the two spherical nuclei are not much frag- 
mented and are found to be more collective than dipole 
and monopole responses. 

Table IIIII gives the theoretical mean energies of 
the GQR compared with values from the systematics 
E(GQR) = 63 A~3 MeV and available measurements 
[36|, 134|. The same energy range, [9-41 MeV], is consid- 
ered to calculate the mean theoretical and experimental 
values [3a, [331 • Contrary to monopole results, theoretical 
GQR energies are close to systematic values but larger 
than experimental ones. In the ^''Mg nucleus, the theo- 
retical energy of the GQR (exhausting 90% of the EWSR) 
is 3.64 MeV larger than the experimental one (exhausting 
only 72 ± 10% of the EWSR) ^. In the ^^Si nucleus, 
the theoretical energy of the GQR (exhausting 91% of 
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FIG. 5: (Color online) Fractions of the quadrupole EWSR in ^^'^^Mg and in ^'^"='°Si. In the deformed nuclei, K'' =0+, 
jK|'^ _ 1+ gjj^ |f(;|T _ 2+ components are indicated in black, green and purple, respectively. Blue curves are obtained from 
folding of QRPA spectra with a Lorentzian distribution. 



TABLE III: theoretical (present work), systematic and exper- 
imental mean energy values of the GQR (in MeV) 



E{GQR) MeV 
Th. Syst. 



Exp. 



^^Mg 


20.62 


22.5 




'^Mg 


20.54 


21.8 


16.9±0.6 


^''Mg 


20.97 


21.3 




2«Mg 


21.41 


20.7 




2«Si 


20.63 


21.3 




2«Si 


20.41 


20.7 


18.54 ± 0.25 


3°Si 


21.44 


20.3 





the EWSR) is 1.87 MeV larger than the the experimental 
one (exhausting only 68 ± 9 % of the EWSR) [H]. The 
mismatch of 10% in ^^Si and 20% in ^"^Mg with experi- 
mental values is probabely the consequence of the value 
of effective mass {m* /m = 0.7) of the DIS interaction 
which is the one giving correct single-particle properties 



in mean-field calculations. Let us note that a reanaly- 
sis of the experimental data in ^®Si nucleus has given a 
higher value of the centroid energy 18.77 ± 0.35 MeV cor- 
responding to 102 ± 11 % of the E2 EWSR in the region 
[9 - 35 MeV] [40|. Keeping in mind that our theoretical 
description does not include any adjustable parameters, 
the agreement with experimental data can be considered 
quite satisfactory. 



C. Octupole response 

Fig. [5] shows the octupole responses in 22-28]y[g g^^^^ 
26-30si. In this figure, K'^ == 0", K'^ = ±1", K'^ = ±2" 
and K'^ = ±3~ components are indicated in black, green, 
purple and red, respectively. As for the quadrupole re- 
sponse, a folding of the QRPA spectra with a 2 MeV 
width Lorentzian distribution has been performed. As 
known for a long time [4l|, the octupole response is 
split into two energy components, the low energy oc- 
tupole resonance (LEOR) generally collective and lit- 
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FIG. 6: (Color online) Fractions of the octupole EWSR in ^^-^^^Mg and ^^'-^''Si. In the deformed nuclei, K'^ = 0~, IK]'' = 1", 
|K|'^ = 2~ and |K|'^ = 3~ components are indicated in black, green, purple and red, respectively. Blue curves are obtained 
from folding of QRPA spectra with a Lorentzian distribution. 



tie fragmented, and the high energy octupole resonance 
(HEOR). In all nuclei studied here, these two energy 
components are well identified, the LEOR ranging from 
5 MeV to 25 MeV and the HEOR from 25 MeV to 
45 MeV. In the oblate ^s-asgj ^nd ^'^Mg, the LEOR is 
typically made of a highly collective [Kj'^ = 3~ state 
around 8 MeV and K-mixed fragments around 15 MeV. 
In the prolate ^^~^^Mg isotopes the LEOR is even more 
fragmented without major contribution from any given 
angular momentum projection K. In deformed nuclei, all 
K-components contribute to the HEOR, which is found 
to be highly fragmented. However, the ordering of the 
different K-components is still found inverted between 
prolate and oblate nuclei: between 30 MeV and 35 MeV 
|K|'^ — 1~ states are those contributing the most to the 
strength in prolate nuclei, whereas in oblate nuclei, the 
strength is due mainly by |K|'^ = 3~ states. Octupole 
responses in the spherical nuclei ^°Si and ^^Mg also show 
LEOR and HEOR components. 



Table IIVI gives the theoretical mean energies of LEOR 
and HEOR. They are calculated in the two intervals 

[5-25 MeV ] and [25-45 MeV ], respectively. Val- 
ues calculated from the energy systematic law of the 
LEOR (E = 30 A-^/^ MeV) are also given. The the- 
oretical values of LEOR are larger than those from sys- 
tematics. They rather follow a "corrected" systematics 
E = 35 A" 1/3 MeV. No systematics exists for HEOR in 
light nuclei. Nonetheless our results are compatible with 
the systematic law E(HEOR) == 110 A ^/^ MeV, obtained 
from experimental data in heavier nuclei [4l| 



VII. CONCLUSION 

A fully consistent microscopic axially-symmetric- 
deformed QRPA approach has been developed and ap- 
plied to light even-even s-d shell nuclei. In these calcula- 
tions, the same Gogny DIS effective force is used in both 
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TABLE IV: QRPA LEOR and HEOR mean energy values 
(in MeV) obtained integrating from 5 to 25 MeV and 25 to 
45 MeV, respectively. Values from the systematics of the 
LEOR are also given. 



(MeV) 


LEOR 




HEOR 




Th. 


Syst. 


Th. 


22Mg 


12.37 


10.70 


34.48 


24 Mg 


13.49 


10.40 


34.83 


26Mg 


12.13 


10.12 


35.72 


2«Mg 


11.18 


9.87 


35.23 


2«Si 


11.88 


10.12 


34.97 


2»Si 


11.54 


9.87 


34.97 


30Si 


11.67 


9.65 


35.63 



HFB and QRPA approaches in all ph, pp and hh chan- 
nels. Results on giant resonances have been obtained in 
22-28]y[g g^^^ 26-30gj igotopes and the influence of the 

ground state intrinsic deformation on strength distribu- 
tions has been discussed. 

In deformed nuclei, theoretical isovector dipole re- 
sponses are found split into two major components, as 
expected. The splitting is found to be correlated with 
specific families of ph transitions. The strength is pre- 
dicted to be almost equally distributed between the two 
energy components. The two groups of transitions corre- 
spond to dipole vibrations along the major and the minor 
axes of the deformed nuclear systems. K = components 



are found in the low energy part of the spectra in the pro- 
late nuclei, and |K| = 1 components in the oblate ones. 

Isoscalar monopole resonances also display a splitting 
in deformed nuclei. Such a splitting is not found when a 
spherical HFB state is used to construct QRPA solutions. 
Isoscalar quadrupole and octupole resonances are found 
to be well fragmented in particular in well-deformed nu- 
clei. Such a fragmentation masks the deformation split- 
ting and the different K broad components are overlap- 
ping. However, LEOR and HEOR components are pre- 
dicted to be well separated, and the ordering of the K 
components as a function of energy is found to be related 
to the sign of the quadrupole moment. Results obtained 
in ^^Mg and in ^*Si show that the Gogny interaction qual- 
itatively reproduces known resonances without resorting 
to any readjustment of parameters. 

These results show that a fully microscopic QRPA ap- 
proach of deformed nuclei is able to provide a wealth of 
structure information and a satisfactory agreement with 
available data. Extensions of the present work to other 
nuclei, in particular exotics ones and to phenomena such 
as pygmy states are under study. 
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